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ABSTRACT 



Currently the only technique sensitive to Earth mass planets around nearby stars (that are too close for microlensing) is the monitoring 
of the transit time variations of the transiting extrasolar planets. We search for additional planets in the systems of the hot Neptune 
GJ436b, and the hot- Jupiter XO-lb, using high cadence observations in the J and Ks bands. New high-precision transit timing 
measurements are reported: GJ436b Tc = 2454238.47898 +0.00046 HJD; XO-lb Tc(A) = 2454218.83331 +0.00114 HJD, Tc(B) = 
2454222.77539+0.00036 HJD, Tc(C) = 2454222.77597+0.00039 HJD, Tc{D) = 2454226.71769+0.00034 HJD, and they were used 
to derive new ephemeris. We also determined depths for these transits. No statistically significant timing deviations were detected. 
We demonstrate that the high cadence ground based near-infrared observations are successful in constraining the mean transit time to 
~ 30 sec, and are a viable alternative to space missions. 

Key words. Stars: planetary systems - Stars: individual: GJ436, XO-1 - Methods: observational 



1. Introduction 

More than 300 extrasolar planets are known up to date. Most of 
them were discovered by radial velocity searches. At this time 
[ nearly sixty transiting systems have been detected. Some of them 
■ reside in crowded fields, making them difi&cult targets for follow 
' up studies. The small amplitude of the transits (<3%) and the 
low probability of the suitable geometric configurations (orbital 
plane along the line of sight, i.e. edge-on orbit) partly explain 
the number of detections. 

Assuming that the host star radius is known, transits allow 
a direct determination of the orbital plane inclination and the 
planet radius. When combined with radial- velocity data, they 
give access to the real planet mass and to the mean planet den- 
' sity revealing its nature: gaseous, icy or rocky. Transits are po- 
tentially sensitive to Earth size planets which cause too small 
movement of the host star to be detected via radial velocity. Of 
course, direct occupations of stars by Earth like planets are also 
difl&cult to detect and obtaining adequate observations is feasible 
only with space missions. 

A numerical investigation of the t ransit timing se nsitivity 
to perturbing plane t s was reported in Agol et al.' (l2 005|) . and in 
iHolman & Murray! (l2005l) . The transit time variations (TTVs) 
depend on the detailed configuration of the system, but some 
fairly robust general predictions are possible for the case of hot 
Jupiters (such as most of the transiting planets known so far): 
(i) TTVs are largest in case of resonant orbits; (ii) TTVs are 
proportional to the mass of the perturbing planet; (iii) TTVs are 
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proportional to the period of the transiting planet; (iv) TTVs vary 
with time (see Eig. 6 in Agol et al. 2005), and in an optimal or- 
bital configuration the interval between two sequential transits 
can change by as much as few minutes. Eor example, an Earth 
mass planet in the system of HD 209458, in a 2: 1 resonance with 
HD 209458b will lead to TTV of order of 3 min. 

The probability for observing long period transiting plan- 
ets is low and all known transiting planets have periods 
<9.2days, with the only exceptions of th e high-eccentricity 
p lanets HD 17156b (Barbieri et al. l2007|), and HD 80606b 
(iMoutou et al.ll2009l) . with orbital periods of ~ 2ld, and ~ 1 1 Id 
respectively. However, systems of multiple extrasolar planets on 
resonant orbits were found by the radial velocity searches (i.e. 
51 Cancri, 3:1 resonance; HD 82943 and Gl 876 - 2:1), showing 
that this is a realistic possibility. The available data show that at 
least 12% of known extrasolar planetary systems have more than 
one planet, and possibly the real fraction is much larger. 

The insuflftcient timing accuracy of individual transits limited 
the planet searches with this method so far: most timing data 
came from small telescopes with large "dead" time between the 
images, in some cases 2/3 or more of the observation duty-cycle. 
Eurthermore, many of these observations require defocusing of 
the telescope (usually because of the large pixel sizes) leading to 
contamination from fainter neighboring star s, sometimes equal 
to the flux of the host star ( Bakos et al.ll2QQ6h . 

The first detailed TTV study (ISteflTen & Agolll2QQ5h used the 
observation of 12 transits of TrES- lb by lAlonso et al.l (l2004 and 
ICharbonne au et al.* (2005) to search for additional planets in this 
system. They found no convincing evidence for a second planet, 
and they can only set an upper mass limit for planets in low or- 
der resonances comparable or lower than the Earth mass, mak- 
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ing these timing data the first that are sensitive to Earth mass 
perturbing planets. 

The possibiHty of having another planet in th e system of 
OGLE-TR- 1 1 3b was studied bv lGillon et all (l2QQ6h - their max- 
imum TTV amplit ude is 43 sec (2.5c r), consistent with < 1-7 
Earth mass planets. iBakos et all (l2QQ6l) derived residuals for in- 
dividual transits for HD 189733b ranging from 0.7 sec to 302 sec 
or 0. 0045-5. 9(T, and the errors of individual transit times are 19- 
150 sec. They are consistent with perturbations from 0.15Mjup 
mass planet at 2: 1 resonance orbit that would remain undetected 
in radial velocity observations. The authors refrain from strong 
statement because th e data are aff'ected by systematic errors. 
Similarly, iDiaz et all 12008) found variations in the period of 
transits of the planet OGLE-TR-111 whose origin has not been 
conclusively determined. 

Space based photometry developed by the MOST team has 
provided accurate values for the transit time s on the HD 189733 
and HD 209458 systems ( Miller-Ricci et al.l>2008a.b) . These re- 
sults have ruled out the presence of super-Earths in the inner 
resonances. Analysis of the transit times for ground-based ob- 
servations of various transiting systems hav e been performed by 
the Transit Light Curve (TLC) team (e.g. Holman et aD 120061: 
IWinn etal.|[2QQ7b . who have obtained accurate timing values 
which have show no strong evidence for the presence of a third 
body in the systems. 

Here we describe the first results from our timing study of 
individual transits with infrared (IR) detectors that allow us to 
obtain imaging with minimum "dead" time for readout (<0.1%). 
By design the IR detectors read out faster than the CCDs be- 
cause the high background forces the usage of short exposures 
and the IR array technology has evolved to achieve reset/readout 
times of order of microseconds, rather than the many seconds 
needed to shift the charges across the CCDs. This gives us the 
following advantages: (i) we observed with unprecedented time 
resolution of -0.1-0.2 sec; the host stars on individual images 
have S/N~50-100, depending on the band and the target bright- 
ness; (ii) we observed bright planet hosting stars without defo- 
cusing (the exception was one transit observation of GJ-436b), 
as it is often the case with the previous studies that attempted to 
use larger telescopes, reducing the contamination from nearby 
sources; (iii) we relied on the ESO timing system that provided 
us with uniform time accurate to better than 0.1 sec, minimizing 
any systematic eff'ects - a crucial advantage over other studies 
that rely on collecting data from various telescopes. Our pre- 
liminary simulations that included only well-behaved Poisson 
photon noise sugg ested a transit timing accuracy of 0.1-1 sec 
(llvanov et al ] |2009I) . 

We apply this high-cadence method to GJ 436b - the first 
transiting hot Neptune planet reported ( Gillon"^^ ^07b), 
hosted by a M2.5V star. The small size of GJ436b leads to a 
challengingly shallow depth of only 0.6%. This planet is partic- 
ularly interesting because it shows large eccentricity which may 
be caused by the gravitational perturbation of a third body in 
the system. The possible presence of a su per-Earth ne a r a 2: 1 
mean motion resonance was proposed by iRibas et aP (I2008D . 
This scenario was recently ruled out, b ut the presence of a third 
body may still be possible (iRibas et al .112 009). We also apply this 
met hod to the Jupiter-mass planet XO-lb, hosted by a Sun-like 
star (iMcCullough et al...2006) . in a ~ 4 d orbit. 

2. Observations and Data Reduction 

The XO-lb data were collected with the SofI (Son of ISAAC) 
instrument at the 3.6-m ESO New Technology Telescope (NTT) 



on La Silla, and with the ISAAC (Infrared Spectrometer And 
Array Camera) instrument at the 8.2-m UTl (Antu) unit of the 
ESO Very Large Telescope on Cerro Paranal, in Visitor Mode, 
and GJ 436b data were collected only with Sofl. All observations 
were carried out in the Fast-Phot cube mode, which produces a 
series of data-cubes with short integration times, and with vir- 
tually zero dead time between integrations, because of the re- 
duced communications between the detector and the instrument 
workstation in this mode. During the observations the detector 
was windowed down to minimize the readout and data transfer 
overheads, with the requirement that the field of view contained 
the target and a reference star of a similar brightness, which 
was used for diff'erential photometry. The windowing allows us 
to reduce the detector integration time to less than 0.01 sec, if 
necessary. The typical overhead is ~ 6 sec per cube of 100- 
2000 frames. A summary of the observing details is presented 
in Tableffl 

2.1. Sofl Observations 

Sofl is the infrare d camera and spectrogra ph at the NTT tele- 
scope on La Silla (iMoorwood et al.|[T998bh . It is equipped with 
a Hawaii HgCdTe array of 1024x1024 pixels, with a gain of 
5.4e"ADU"^ and a readout noise of 2.1 ADU. Its detector 
shows a non-linearity of less than 1.5% below 10,000 ADU, in a 
Correlated Double Sampling readout. The imaging mode has a 
pixel scale of 0.288 arcsecpix"^ 

The observations of GJ436b were carried out the night of 
May 17, 2007, in poor weather conditions. We applied a lin- 
earity correction to the data based on calibration data obtained 
on May 14, 2007, where a 4-order polynomial was fitted to the 
deviation from linear detector responsq^- Despite of defocusing 
the telescope in this run, some images show pixels values above 
the correctable 1 8 000 counts level, and they were omitted from 
the resulting light curve, such that the final sample spans 35 266 
points, covering -146 minutes. 

We observed two transits of XO-lb with Sofl, during the 
nights of April 27 (run A), and May 1 (run B), 2007, both in 
the / band, covering 228 min and 303 min, respectively. 

2.2. ISAAC Observations 

ISAAC is an infrared camera and spectrograph locat ed at 
the Nasmyth B focus of UTl (iMoorwood et al.1 Il998ah . For 
our observations we used the long-wavelength arm, which 
equipped a 1024x1024 Aladdin array, with a pixel scale of 
0.148arcsecpix"\ a gain of 8.7e"ADU"\ and a read noise 
of 4.6 ADU. The readout mode was Double Correlated Read 
Low Bias. This detector is linear at 90% for signal below 
16,000 ADU. 

We observed two transits of XO-lb, during the nights of May 
1 (run C), 2007, in the Ks band, and May 5 (run D), 2007, in the 
/ band, using the ISAAC /-hBlock filtefl. These observations 
cover 277 min and 257 min, respectively. 



^ More information can be found at the Sofl web page 
http://www.eso.org/sci/facilities/lasilla/instrunients/sofi/inst/Linearity.htmTl 

^ The J filter suff'ers from a red leak. Normally, it is eliminated by the 
sensitivity cut-off" of the Hawaii detectors at ~ 2.5/um, but the Aladdin 
detector used in the long-wavelength arm of ISAAC is sensitive to 
~ Sjim. A blocking filter is added to eliminate the leak. The overall 
transmission of the /+Block filter is similar to that of the "standard" / 
filter. 
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Target 


Run 


Date 


Instrument Window size (px) 


DIT (sec) 


N. Frames 


Filter 


Apertures (px) Inner annuli (px) 


GJ436b 


A 


May 17, 2007 


SofI 


872 X 132 


0.239 


36058 


Ks 


23 


35 


XO-lb 


A 


April 27, 2007 


SofI 


352 X 152 


0.8 


17 066 


J 


11 


19 


XO-lb 


B 


May 1, 2007 


SofI 


892 X 180 


0.8 


20100 


J 


8 


20 


XO-lb 


C 


May 1, 2007 


ISAAC 


64 X 544 


0.17-0.13^ 


110161 


Ks 


7 


10 


XO-lb 


D 


May 5, 2007 


ISAAC 


64 X 544 


0.08 


186095 


/-hBlock 


21 


21 



^ The DIT (Detector Integration Time) was reduced during the observations to avoid reaching the non-linear limit of the detector in improving 
sky conditions. 



2.3. Data Reduction 



Table 2. Transit light curves. 



Standard infrared data reduction steps were applied: flat field- 
ing, and dark subtraction. However, since the observations were 
obtained in stare mode, i.e. with no jittering, the sky was not sub- 
tracted with the usual method used for infrared data. Instead, we 
estimated the sky level by measuring the flux in circular annuli 
centered on the target and the reference star. Note that PSF fit- 
ting was not possible because of the defocusing of the GJ436b 
observations and even if it would not have been the case, the 
targets are the brightest sources in the field, and the seeing varia- 
tion did not allow to create a PSF model from the previous and/or 
next frames. The fundamental limit of how bright the reference 
star could be comes from the maximum size of the detector win- 
dow. We always select as a reference source the brightest avail- 
able star in the field of view, and if it is fainter than the tar- 
get, it dominates the noise of the final light curve, as in the case 
of GJ436b. The data reduction was carried out with the IRAlEl 
package DAOPHOT. 

The final light curves were divided by a linear polynomial of 
the form Correction = a-\-bxT, where T is the time, calculated 
with the out-of-transit points of each light curve, to normalize 
the light curve, and correct any smooth trend due to atmospheric 
variations. 

The five light curves obtained are presented in Table|2l where 
we only show a small fraction of the complete light curves as a 
format example. The complete light curves are available in the 
electronic version of the Journal. 



3. Analysis 

We focus on measuring the central transit time (Tc) and depth 
(d) . For GJ 436b we ado pted the stellar parameters determined 
bv iGiUon et all (l2007bl) : Te/f = 3500 K, logg = 4.5, and 
[Fe/ H1 = Q-Q- The rest of the system parameters were taken 



HJD 



Flux 



Error 



from 
from 



TbrresI (12007^ XO-l b, we adopted stellar parameters 
McCullou ghetaD (12006 ): Te/f = 5750 K, logg = 4.53, 
and [ M/H] = 0.058. Planetary and orbital parameters were taken 
from iHolman et al.l ^2006). The transit length was calculated 
from the known orbital parameters and the depth was measured 
as ratio of the out-of-transit to the in-transit flux. The next step 
was to create a light curve m odel according to the prescrip- 
tion of Mandel & A20I (2002), assuming throughout the paper 
a quadratic limb-darkening law , where the limb-darkening co- 
efficients were taken from IClaretl (|2QQQ|) for the adopted stellar 
parameters, and the used band passes. 



^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



GJ 436b 
2454238.456329 0.986084 
2454238.456332 
2454238.456335 
2454238.456337 
2454238.456340 



0.995099 
1.005212 
1.018050 
1.015619 



xo- 

2454218.911636 
2454218.911645 
2454218.911654 
2454218.911663 
2454218.911673 



lb Run A 
1.003098 
1.017395 
1.017412 
0.991108 
0.988897 



0.005644 
0.005565 
0.005784 
0.005828 
0.005785 



0.003059 
0.003063 
0.003075 
0.002969 
0.003064 



The final step was to fit this light curve to the observations 
minimizing the statistics: 



X 



fi-f 



(1) 



where f^^ is the observed flux with an uncertainty o-f, and fi 
is the expected flux obtained from the model. Here, the central 
transit tim e (7^) was th e only free parameter . The minim ization 
method of lBrentI (Il973l) . as implemented in lPress et al.1 (I1992D . 
was used. 



3.1. G J 436b 

The transit depth we measured for GJ 436b is d = 0.64 ± 0.03%. 
This corresponds to a planet- to- star size ratio: p = 0.0 8 2 + 0.0 02 
which is in good agreement with the iGiflon et al.l (l2007allbl: 
p = 0.0829 ± 0.0 043 and p = 0.082 ± 0.005, respectively), 
and marginally with lDeming et al.l ([2007; p = 0.0839 ± 0.0005). 
For the values given above, we adopted the corresponding limb- 
darkening coefficients in ^: a = -0.0677, b = 0.3665. 

We fitted the model to the observed light curve to obtain the 
transit midpoint Tc = 2454238.47898 ± 0.00046 HJD (~ 39 sec 
error). The observed light curve for GJ436b is shown in Fig.[T] 
(top panel), together with a light curve binned to 30 sec time res- 
olution, and the fitting model (middle panel), and the residuals 
of the binned curve (bottom panel). The binned light curve was 
used in this analysis. The time and the flux values of a bin are 



^ We assume the same limb darkening coefficients for K and Ks be- 
cause the effective wavelengths of the two filters are similar. 
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0.92 

1.010 
^ 1.005 
^ 1.000 
I 0.995 
^ 0.990 
0.010 
0.005 
0.000 
-0.005 
-0.010 



SofI 

ay.^ = 15.0 mmag 
aJ = 35 sec 




H 1 1 \ \ 1 1 1 1 \ 1 1 1 1 \ \ 1 1 1 1 \ \ 



Oy. = 1.4 mmag 




H 1 1 \ \ 1 1 1 1 \ 1 1 1 1 \ \ 1 1 1 1 \ \ 




238.46 238.48 238.50 238.52 238.54 238.56 

HJD - 2454000 

Fig. 1. Photometry of the GJ436b transit obtained with Sofl. 
Top panel: the normalized light curve, with a integration time 
of 0.239 sec. The central time of the transit is marked at HJD 
2454238.47898. Error-bars are omitted for the sake of clarity. 
The out-of-transit r.m.s. flux on the unbinned curve (crKg), and 
the timing error (0-^^) are also shown. Middle panel: the binned 
curve, and the best fitting model. The bin width is 30 sec. Error 
bars represent the flux r.m.s. within each bin. Note the difl'erent 
Y-axis range on the panels. The out-of-transit r.m.s. flux (ctks ) of 
the binned curve is also presented. Bottom panel: The residuals 
of the best fitting model to the binned curve. 



the average of the times and fluxes of all measurements within 
the bin, and the flux error is their r.m.s. 

We performed a bootstrapping simulation to calculate the 
time of transit uncertainty. The set of residuals of the best-fitting 
model was shifted a random number of points in a circular way, 
and then added to the model light curve, constructing a simu- 
lated light curve with the same point-to-point correlation as the 
observed light curve. This procedure takes into account the cor- 
related noise in our analysis. Then, we calculated the center-of- 
transit for the new curve, as described above. This procedure 
was repeated 10000 times, and the l-cr width of the resulting 
distribution of timing measurements was adopted as the error of 
the timing. The l-cr error value is weakly dependent on binning. 
Here we choose a bin size of 30 sec because it is a good compro- 
mise between the number of points included in each bin and the 
final number of points in the light curve. 

Many follow-up observations of GJ 43 6b has been carried 
out with both Spitzer Space Te lescope jG illon et al. 2007al: 
iDeming et af][2Q07: Dem orv et al.l l2007l: [Southworth 20081) and 
the Hubble Space Telescope (Bea n & SeifahrtI l2008l) . and re- 
cently ground based observations have given a timing precision 
comparable to space based observations ( Alonso et al. 2008; 
fShporer et al.ll2008l) . They are all Hsted in Table[3l 

Considering the newest literature data, and our measurement 
we recalculated the ephemeris of GJ436b by fitting a weighted 
linear relation, to obtain a period P = 2.6438986 ± 0.0000016 d, 
and a "zero transit" epoch Tc(E = 0) = 2454222.61588 ± 
0.00012 HJD. The new epoch presented in this work is in ex- 
cellent agreement with these ephemeris, and they in turn agree 



40 60 
Epoch (HJD 



80 

2454222) 



Fig. 2. The observed minus calculated time-of-transit diagram 
for the diff'erent timing values in the literature, as a func- 
tion of the observed epoch, for GJ436b. Diff'erent symbols 
repre sent difl' erent literature sources. The open squares are 
from IShporer e t al. (2008). The open tri angle is an average of 
the value fr o mlGillon et al.1 (|2007a), De ming et al.1 (l2007hT and 
ISo uthworthI (l2008l) . with the corr ections given by iBean et al.l 
( 2008|). The solid t riangles are fromlBean & SeifahrtI (|2008|). The 
open circle is from Alo nso et al. 1('2008). The solid squares come 
from Ribas et al. ( 2009^). Finally, our timing estimate is drawn 
with a solid circle aiE = 6. All measurements are listed in Table 

El 



Table 3. GJ436b timing measurements. 



Tc (HJD) 


o-Tc (d) 


Epoch 


Reference 


2454222.61612 


0.00037 





(1) 


2454225.26029 


0.00038 


1 


mi 


2454238.47898 


0.00040 


6 


This work 


2454243.76657 


0.00026 


8 


(1) 


2454246.40982 


0.00040 


9 


(1) 


2454251.69956 


0.00030 


11 


(1) 


2454280.78165 


0.00017 


22 


(2,3,4)^ 


2454439.41607 


0.00068 


82 


(5) 


2454444.70385 


0.00093 


84 


(5) 


2454447.34757 


0.00080 


85 


(5) 


2454463.20994 


0.00089 


91 


(5) 


2454468.49911 


0.00094 


93 


(5) 


2454505.51379 


0.00050 


107 


(6) 


2454534.59584 


0.00015 


118 


(7) 


2454558.39010 


0.00063 


127 


(6) 


2454587.47447 


0.00061 


138 


(6) 


References: (1) - Shporer et al. ('2008D 
(2007a); (3) - Demins et al. (20oC 
(2008); (5) - Bean & Seifahrt (120081) 


; (2) -IGillon et al, 
(4) -ISouthworth 
(6) -iRibasetaL 



(2009); (7) - Alonso et al. (2008 

^ Average of the two values given for this epoch. 

^ Average of the three values given by different authors 
for this epoch, considering the correction given by Bean 
et al. (2008) 



with those of 'Ribas etall (|2009l), and iBean & Se ifahrtI (l2008b. 
Figure |2] shows the residuals of the fit of the new ephemeris as a 
function of the observed epoch for the available timing values in 
the literature, and our timing value at epoch E = 6. 

The data show some TTVs of up to 98 sec. - smaller than the 
predicted deviations of order of a few min utes for a 1-10 Eart h 
mass companion on a resonant 2:1 orbit (lAlonso et al.ll2008l) . 
However, this deviations are consistent with zero, within their 
respective uncertainties. Further observations with higher accu- 
racy are necessary to constrain better the properties of this sys- 
tem and to address the question if it bears other planets. 
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Table 4. XO-lb timing measurements. 



Tc (HJD) 


(Ttc (d) 


Epoch 


Reference 


2453127.03850 


0.00580 


-173 


(1) 


2453142.78180 


0.02180 


-169 


(1) 


2453150.68550 


0.01060 


-167 


(1) 


2453154.62500 


0.00260 


-166 


(1) 


2453158.56630 


0.00340 


-165 


(1) 


2453162.51370 


0.00250 


-164 


(1) 


2453166.45050 


0.00250 


-163 


(1) 


2453170.39170 


0.00370 


-162 


(1) 


2453229.51430 


0.00450 


-147 


(1) 


2453237.40430 


0.00320 


-145 


(1) 


2453241.34100 


0.00670 


-144 


(1) 


2453808.91700 


0.00110 





(2) 


2453875.92305 


0.00036 


17 


(3) 


2453879.86400 


0.00110 


18 


(3) 


2453883.80565 


0.00019 


19 


(3) 


2453887.74679 


0.00016 


20 


(3) 


2454218.83331 


0.00114 


104 


This work (A) 


2454222.77539 


0.00036 


105 


This work (B) 


2454222.77597 


0.00039 


105 


This work (C) 


2454226.71769 


0.00037 


106 


This work (D) 



References: (1) - Wilson 

iMcCullough ^tall (12006) : (3) 



et al. (2006); (2) 
Holman et al.. (i2006) . 



3.2. X0-1b 

The XO-lb light curves were analyzed in a similar way as for 
GJ436b. The quadratic stellar limb-darkening coefficients uti- 
lized here were: a = 0.00592 and b = 0.34954 for the Ks light 
curve, and a = 0.10923 and b = 0.35938 for the /-hBlock light 
curve. 

After fixing the system parameters we calculated the transit 
midpoint for the 4 light curves separately. The errors were cal- 
culated with the same bootstrapping technique described above. 
The transit timings calculated here are shown in Tabled] 

Interestingly, the uncertainties of the runs A and B decrease 
significantly, if the;^^ is calculated only over the ingress and the 
egress phases: 0.00012 d and 0.00009 d (10 and 8 sec), respec- 
tively. The bootstrapping simulation was also done over these 
parts of the light curve. Apparently, using only the ingress and 
egress excludes some of the systematic eff'ects that occurred dur- 
ing the rest of the transit, and were reflected in the error distri- 
bution. Therefore, we consider the errors given in Table H] to be 
upper limits of the uncertainties. The errors of the Runs C and 
D remained virtually unchanged: 0.00044, in both cases, most 
likely because they were obtained in poor weather conditions 
and the errors are dominated by the reduction of flux during the 
periods of poor atmospheric transmission. 

The final light curves are shown in Fig. [3j with a 20 sec bin 
width version to easily show the best fitting model. Run A shows 
some systematics that could not be corrected, so this curve was 
only used to get timing value, and not planetary parameters. 

We fitted a weighted linear relation to the timing values listed 
in Table m to get the predicted ephemeris for the transits of XO- 
lb. This new ephen ieris correct the lon g term diff'erence in the 
ephemeris given by iMcCullough et aP (2006) and Wilson 'etaD 
(120061) . Our fit gives us: P = 3.9415128 ± 0.0000028 d, and a 
"zero transit" epoch Tc(E = 0) = 2453808.91682 ± 0.00013 
HJD. In this calculation, we use the weighted average of runs 
B, and C, which spans the same epoch. This calculations are 
shown in Fig. IH Note that runs C and D, carried out with the 
bigger telescope, but under worse weather conditions yield less 
accurate timing measurements than transits observed with the 
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Fig. 4. The observed minus calculated time-of-transit diagram 
for the diff'erent timing values in the literature, as a function of 
the observed epoch for XO-lb. For clarity the diagram was di- 
vided in two diff'erent scale panels (minutes in the left panel, 
and seconds in the right panel). Diff'erent symbols r epresents dif- 
ferent literature sources. The pentagons are from [ Wilson et"aD 



(l2QQ6h . The squares are the values from iHolman et al 



The open triangle is the value from IMcCullough et al 



Our timing values are shown with solid triangles. All measure- 
ments are listed in Table (H 



smaller telescope but under better weather conditions (i.e. run 
B), demonstrating the impact of the weather on the high preci- 
sion photometry required to detect the transiting planets. 

The amplitude of the resulting TTVs, show no evidence for a 
perturbin g third b ody in the system, in agreement with the results 
of .Holman et al.1 (120061 



4. Conclusions 

Here, we present new ground based high cadence near-infrared 
observations of one transit of the hot-Neptune GJ 436b, and four 
transits, spanning three epochs, of the hot- Jupiter XO-lb. 

We achieve transiting timing accuracies of about 30 sec for 
individual transits. The uncertainty is dominated by systematic 
eff'ects, and greatly exceeds the few second errors predicted by 
photon noise dominated observations. We find no significant ev- 
idence for perturbations of the orbital motion of GJ436b nor 
XO-lb by other bodies in the system. Of course, a proper test of 
this hypothesis will require monitoring of multiple transits with 
the same or even higher accuracy. 

We demonstrate that the ground-based high-cadence obser- 
vations of transiting extrasolar planets is an excellent technique 
for constraining the parameters of extrasolar planetary systems 
because of the statistical significance of the obtained timing 
measurements. The timing precision is comparable with the 
space-based observations, making this method a good alternative 
to the space mission with their high cost and limited life-time. 
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Fig. 3. Photometry of the XO-lb transits obtained with SofI and ISAAC. Each panel shows (from top to bottom) the original light 
curve with a 10 sec bin size version with the best-fitting model, and the residuals of the fit. Top-left panel: Run A, observed with 
SofI, in the / band. Top-right panel: Run B, observed with SofI, in the / band. Bottom-left panel: Run C, observed with ISAAC, 
in the Ks band. Bottom-right panel: Run D, observed with ISAAC, in the /-hBlock band. Each panel gives the out-of-transit r.m.s. 
flux of the unbinned curve (upper sub-panel), and of the binned curve (middle sub-panel). The timing error (cr^^) for each transit 
is printed. Note that for runs A and B we also give, in brackets, the timing error obtained if only the ingress-egress phases are 
considered. See more details in Sect. 13.21 



